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Abstract. Some types of coupled Korteweg de-Vries (KdV) equations are derived from 
an atmospheric dynamical system. In the derivation procedure, an unreasonable y- 
average trick (which is usually adopted in literature) is removed. The derived models 
are classified via Painleve test. Three types of r-function solutions and multiple soliton 
solutions of the models are explicitly given by means of the exact solutions of the usual 
KdV equation. It is also interesting that for a non-Painleve integrable coupled KdV 
system there may be multiple soliton solutions. 



1. Introduction 

The single component KdV equation has been widely used in various natural science 
fields especially in almost all branches of physics. For instance jT], the KdV equation 
describes, in a general form, competition between weak nonlinearity and weak disper- 
sion, while the nonlinear Schrodinger (NLS) equation describes the same competition for 
envelope waves (see, for example, the Introduction in [2]). Some other integrable equa- 
tions such as the sine-Gordon (SG) equation, the Kadomtsev-Petviashvily equation, the 
so-called three- and four- wave systems are universal as well. 

Some kinds of the coupled KdV equations had also been introduced in literature such as 

one describing two resonantly interacting normal modes of internal-gravity-wave motion 

in a shallow stratified liquid In principle, many of other coupled KdV equations are 

introduced mathematically because of their integrability pQ. 
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In section II of this paper, we derive some new types of coupled KdV equation systems 
with some arbitrary parameters from a two-layer fluid model which is used to describe 
the atmospheric phenomena by using long wave approximation. 

Whence the coupled KdV systems are derived, an important problem is how to solve 
them. To get more exact solutions, one hopes to figure out the integrable ones. In section 
III of the paper, we use the well known Painleve test classification to find out the Painleve 
integrable ones for special types of selections of the parameters. 

For some special types of coupled KdV systems, one can find some types of exact 
solutions by modifying the solutions of the usual KdV equation. Some concrete solution 
examples, especially, the r-function solutions and the multiple soliton solutions are given 
in section V. The last section contains a short summary and discussion. 

2. Coupled KdV equations derived from atmospheric dynamics 

It is known that various integrable models can be derived from fluid dynamics. Simi- 
larly, in this section we use a two layer fluid model, 

qu + J{ipi, + = 0, (1) 

q2t + J{^2, qz} + Pihx = 0, (2) 

where 

qi = i>ixx + i>iyy + F(il)2-i)i), (3) 

q2 = ^2xx + ^2yy + - 1p 2 ), (4) 

and J{a, b} = a x b y — b x a y , as a starting point to derive two component KdV equations 
by using the multiple scale approach and the long wave approximation. 

In (JIJ)-(|IJ), F is the weak coupling constant among two layer fluids and (3 = j3o(L 2 /U), 
(3q = (2a; /ao) cos0 o , in which a is the earth's radius, u is the angular frequency of the 
earth's rotation and 0o is the latitude, U is the characteristic velocity scales. The deriva- 
tion of the dimensionless equation (0) and (J2J) is based on the characteristic horizontal 
length scale L = 10 6 m and the characteristic horizontal velocity scale U = 10 _1 m/s 0. 
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More especially, if we take /3 — 0, the system and is reduced to the usual coupled 
Euler equation which is suitable to be used to describe two-layer inviscid fluid. In other 
words, all the results obtained in this paper are valid for general two-layer inviscid fluid. 

Under the long wave approximation in the x direction, in order to derive the KdV type 
equations, the stream functions ip\ and ip2 should have the form 

A = <frio(y) + <f>ii e ( x ~ c ot), V, e 3 t) = iO + <j>i(X, y, T) = (f) i0 + (pi, i = 1, 2, (5) 

where e is a small parameter. It is reasonably considered that the parameters F and /3 
are in order e and e 2 respectively that means the coupling among two layers is week and 
the effect of the rotation of the earth is much smaller. Thus we set 

F = F e, (3 = (3 e 2 . (6) 

Now, we expand the shift stream functions 0,, i = 1, 2 as 

0! = e<f> n {X, y, T) + e 2 12 (X, y, T) + e 3 13 (X, y, T) + 0(e 4 ), (7) 
2 = e0 21 (X, y, T) + e 2 22 (X, y, T) + e 3 23 (X, y, T) + 0(e 4 ). (8) 

Substituting ©-© to JU)-© yields 

[(010?/ - C)d yy - <j)\Qyyy\ <\>\\X^ + {[(010?/ ~ c)3yy ~ filQyyy] 012X + F ((f) 10 y - C)0 2 1X 
+ [F (C - (f)20y) + 4>nyyy}4>UX - <j>Uy<j>llyyx} ^ + {[(010?/ ~ c)d yyX ~ 4>10yyydx} 013 
— 012jy011j,j/X — 4>lly4>V2yyX + (010?/ — C O )(-F o 22 + (pnxx)x + 4>UyyT ~ -^0021X011?/ 
+ [ ( f>12yyy + -^0021?/ + WllX + [F (co - fooy) + <Pllyyy]<Pl2x} ^ + 0(e 5 ) = 0, (9) 

and 

[(0202/ ~ C )dyy ~ V 0yyy] 021XC 2 + {[(020?/ ~ c)d yy - 020?/?/?/] 022X + -Fo(02Oy ~ c)011X 
+ [F (C Q - 010?/) + 021?/?/?/]021X - 021?/021?/?/x} ^ + {[(020?/ ~ c)d yyX ~ <j>20yyydx} 013 
— 022?/021?/?/X — 021?/022?/?/X + (020?/ — Co)(-F O 0i 2 + 021Xx)x + 021?/?/T ~ -?O011X021?/ 

+ [022?/?/?/ + -^0011?/ + /?o]02ix + [Fo(tk) - 010?/) + 02i?/i/?/]022x} e 4 + 0(e 5 ) = 0. (10) 
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Vanishing the e 2 terms of © and ([1U|). we have a special solution 

11 = A 1 (X, T)Bi(y) = A 1 B 1 , (11) 
21 =A 2 (X, T)B 2 (y) = A 2 B 2 , (12) 

with i?! and B 2 linked to 0i O and 2 o by 

tW-Si - B ly <f) Wy + d = 0, 0io = /7 + c y, (13) 

and 

V yyB 2 - B 2y (p 20y + C 2 = 0, 020 = V + c y , (14) 

respectively with arbitrary constants Ci and C 2 . 

By using the relations ([TT]) - ()14j) , the equations obtained by vanishing the e 3 orders of 
(JHJ) and fjlOj) and integrating once with respect to X become 

2<f) l0y (Bxdyy - B lyy ) 12 + B 1 [b n Al - 2F {B 1 <p 20y A 1 + B 2 <p l0y A 2 )} = 0, (15) 
2<f) 20y [B 2 d m - B 2yy ) 22 + B 2 [b 2l A\ - 2F {B 2 <p l0 yA 2 - B^Ax)] = 0, (16) 

where the integrating functions have been dropped away and 

b\\ = B\B\ yyy — Bi y Bi yy , b 2 i = B 2 B 2yyy — B 2y B 2yy . (17) 

It is readily verified that 

12 = (B 3 Al + B 4 A 2 + B A 1 )B ll 22 = (B 5 A 2 2 + B 6 A, + B 7 A 2 )B 2 , (18) 

with i?o, -B3, -B4, i?5, i?6 and .87 being functions of y and determined by 
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solves the third order equations (fTHj) and fjTHfl . 



Because of (fTTj) . (fT2Jl and IJIHjl . the fourth order of © and (HUJ) become 

/i ipyy ~ &i l Bi yy ) 0i3X + BiyyAixT + fiBiAixxx + F Q (giBiB 4: — fiB 2 B 7 )A 2X 
+2f 1 F B 2 B 5 A 2 A 2X - (A)5i - F o9l B B 1 + F f 1 B 2 B 6 )A lx + B i b u (A 1 A 2 ) x 



A X A 



IX 



and 



+ 2b n B — 2F g 1 BiB 3 H — 091 ( 1 2 — - — B 2 B ly + B x B 2y 

J\B 2 V /i #i 

+ 7^2 [ & n(65 3 /i + 3/iS ltf - Cl ) - /iSifeuJ A\A 1X = 0, (23) 



9i (9 yy — B 2 l B 2yy ) 023X + B 2yy A 2X T + 5 , i-B 2 A 2 xxx + F (giB 1 B — fiB 2 B 6 )A lx 
+2g 1 F B 1 B 3 A 1 A lx - {f3 Q B 2 + F f 1 B 7 B 2 - F g l B 1 B 4 )A 2X + 5 6 6 21 (^iA)x 



2b 2 \B 7 — 2F fiB 2 B 5 H — f — — - — 1 2 — B\B 2y + B 2 B\ y 

9i Bi V £i A 

+—j [6 2 i + 3^1^ - dO - <7i£2&2i y ] A^sx = 0. (24) 
^#1 

In the usual studies to solve (j2HJ) and P^J) type equations, especially in the atmospheric 
and ocean dynamics, one would take 0i3 and 023 as zero. However, whence 0i3 and 023 
are taken as zero, there may be a non-consistence problem because Ai and A 2 are only the 
functions of X and T while the coefficients of f!23|) and (}2*4"j) are explicitly ^/-dependent. 
In generally, equations (|2l3|) and are not consistent except that all the ^/-dependent 
coefficients are proportional each other up to constant level. The detailed analysis of 
the equations (j2^|) and (j2H) with 13 = 2 3 = tells that it is impossible to select the 
functions Bo, B\, B 7 , Uq and Vq which are proportional each other and satisfy the 
equations ((EH), (JHJ), (flH|) - (f2~T]) at the same time. To avoid this kind of inconsistency in 
the traditional literature, an unreasonable and unclear procedure is usually made, taking 
a y average by integrating the inconsistent equations with respect to the variable y from 
Vi to y 2 . 

Nevertheless, it is possible to get some consistent and significant solutions from ()23|) 
and (|2^j) by taking nonzero i3 and 2 3- In this paper, we only give out a possible selection 



of and 023 to derive coupled KdV type equations for the quantities A\ and A 2 . 
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It is straightforward to verify that if we take 

013 = n J A lx A 2 dX + r 2 A\ + r 3 A\ + r±A x + r^A x A 2 + r e A 2 2 + r 7 A 2 + r 8 A lxx , (25) 

<P23 — Si 

/ A lx A 2 dX + s 2 A\ + s 3 A\ + s 4 A 2 + s 5 A t A 2 + sqA\ + s 7 A 1 + s 8 A 2X x, (26) 

with 

s l = B 2 f f SMWdy", i = l, 2, .., 8, 

i2i = - aiBlBl ™ , ^2 = ^k[Bihb lly + 6 n (ci - 3/iS lw - 6S 3 /f)], 

Jl OJl 

R > - ~vT [ W3 + — ftA 7? j - 7T ls/!/ + buBo) + HET' 

B B 
Ra = — t-(PqBi — FqBqBiQx + F f 1 B 2 B 6 ), R 5 = — — (a 3 B lyy + b n B 4 ), 
Ji h 

Re = -^[(«2 - a 5 )B lyy + F hB 2 B b \, R 7 = ^(g.B.B, - f\B 2 B 7 ), 
Ji Ji 

B B 

Rs = --r( a 4Bi yy + fiBx), S 8 = (8 4 B 2yy + giB 2 ), 

Ji 9i 

Si = 5lB2B2y \ S 2 = ^[B 29l b 21y + b 21 {d, - 3 9l B 2y - 6B 5 gf)], 
9i Q9i 

BjFpfi f B x B 2y — B 2 Bi y dA CxFpBj B 2 

S4 = (PoB 2 — F B 7 B 2 f 1 + F g 1 B 1 B 4 ), S 5 = — (5 3 B 2yy — b 21 B 6 ), 

9i 9i 

S 6 = —[(6 2 - 5 5 )B lyy - F o9l B 1 B 3 ], S 7 = ^(g^Bo - f\B 2 B Q ) 
9i 9i 

for arbitrary B\ and B 2 , then Ai and A 2 satisfy the following coupled KdV system 

A 1T + a t A 2 A lx + {a 2 A\ + a 3 A x A 2 + «4^ixx + a 5 Al) x = 0, (27) 
A 2T + M 2 ^ix + {5 2 A\ + 5 3 A ± A 2 + 5 4 A 2XX + 5 5 A 2 2 ) X = 0, (28) 

where ten constants {a iy 5 iy i — 1, 2, 3, 4, 5} are arbitrary. 

Now the important question is how to get some types of exact solutions of the coupled 

KdV system. Before giving out some concrete solutions, we try to make a Painleve 
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classification at first. That means we are going to give some constraints on the parameters 
{aj, Si, i = 1, 2, 3, 4, 5} such that the solutions of the model are single valued with 
respect to an arbitrary singular manifold. 

3. Painleve classification of the coupled KdV system 

The Painleve test is one of the best way to study nonlinear systems. In this section, we 
take a standard Painleve test by using the Kruskal's simplification for the coupled KdV 
system. 

To pass the Painleve test, four steps are required. The leading order analysis, the 
resonances determination, the test of the primary branch and the test of the secondary 
branches. 

The leading order analysis for the coupled KdV system ()27|) and (f^Hj) around the arbi- 
trary manifold shows us that there are two possible cases: 
Case 1. 

Ai ~ — A 2 ~ ^ (29) 
In this case the parameters {aci, 5i} and {uo, vo} are linked by 

2a 5 ul + 2a 2 vl + (2a 3 + ai)u t>o + 12a 4 u = 0, 

25 5 v% + 25 2 Wo + (25 3 + 5i)uqv q + 125^ = 0. (30) 

Case 2. 

Ai ~ — , A 2 ~ — (31) 

or equivalently 

A U ° A V ° 

which will be not considered since the exchange symmetry {Ax, A 2 , ai, 5i} <-» {A 2 , A\, 
5i, for the coupled KdV system (|2*7j) and (J2EJ). 
The case (13 lj) appears only for 

5 2 = 0, 5 4 = — (25! + 35 3 ), u = -6—. (32) 
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The resonance analysis for the first case shows us that the resonant points are 
located at 

- 1, 4, 6, ji, j 2 , ja = 9 - ji - j 2 , (33) 
where ji, j 2 and j 3 are three roots of 

u v 8 4 a 4 (j - 9)j 2 + [(U5iV u - u 2 (v 8 1 + 8 3 v + 28 2 u ))a 4 - v%8 4 (2a 2 v + u a 3 )]j 

+ (245 4 w Mo + 2ul(A8 2 u + v 8i + 28 3 v ))a 4: + 2fo^ 4 (M ai + 2u a 3 + ^wo) = (34) 

for the variable j. Apart from the equivalent decoupled case for both A\ and A 2 satisfy 
completely decoupled KdV equations, the positive integer conditions for the resonant 
points lead to the following only nonequivalent subcases (i) j\ = j 2 = 0, j 3 = 9, (ii)ji = 
0, h = 1, 33 = 8, (iii)ji = 0, j 2 = 2, j 3 = 7, (iv)jx = 0, j 2 = 3, j 3 = 6, (v)ji = 0, j 2 = 
4, js = 5, (vi)ji = j 2 = 1, j 3 = 7, (vii)ji = 1, j 2 = 2, j 3 = 6, (viii)ji = 1, j 2 = 3, j 3 = 5, 
(ix)ji = 32 = 2, j 3 = 5, (x)j 1 = 2, j 2 = 3, j 3 = 4. 
The resonance analysis for the second case (}3T|) gives that the resonances will appear 

at 

-1, 0, 4, 6, jt, 32 = 6-3u (35) 

where j\ and j 2 are two solutions of 

j(28 1 + 38 3 )(j - 6) + 275 3 + 225 x = (36) 

for the variable j. It is clear that the positive integer conditions for the resonance points 
lead to the following four nonequivalent subcases (a) j\ = 0, j 2 = 6, (b)ji = 1, j 2 = 5, 
(c)ji = 2, j 2 = 4, and (d)ji = 3, j 2 = 3. 

To check all the resonances for the subcases (i)-(x) and (a)-(d) yields the possible 
Painleve integrable models under some conditions for the model parameters and 5j. 
Here we just list the final nonequivalent models but omit the detail and tedious analysis 
because the procedures are standard. The result shows us that there are only six Painleve 
integrable subcases of the coupled KdV system (|27j) and (|2*Hj) . 
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P-integrable model 1. 



Ait + (a 1xx + ^(c 2 - ci - c lC2 )^ + c 1 A 1 A 2 - ^A 2 ^j = 0, 

A?r+ (a 2X x + ^(ci -c 2 - + - ~cic 2 A^ =0, (37) 

where ci and c 2 are arbitrary constants. For this type of coupled KdV system (|3*7|). there 
are three branches with the resonances { — 1, 2, 3, 4, 4, 6}, { — 1, 2, 3, 4, 4, 6} and 
{ — 1, —1, 4, 4, 6, 6} respectively and all the resonance conditions satisfied. 
P-integrable model 2. 

Ayr + [A lX x - (c + 3)(c + 6) A* - c 2 A 2 2 } x + 2c[(c + 6)A 1X A 2 + (c + 3)AxA 2 x] = 0, 
A 2T + [^2xx - c(c - 3)A 2 2 - (c + 3) 2 Al} x + 2(c + 3) [cA 2 A lx + (c - 3) A, A 2X ] = 0, (38) 

where c is an arbitrary constant. For the model system there is only one branch with 
the resonances located at { — 1, 1, 2, 4, 6, 6} and all the resonance conditions satisfied 
identically. 

P-integrable model 3. 

A 1T + (A 1XX +Al + A X A 2 ) X = 0, 

A 2T + {A 2XX + A 2 2 + A X A 2 ) X = 0. (39) 

In this case, the resonance points are {—1, 0, 4, 4, 5, 6}. 
P-integrable model 4. 

A lT + [Aixx + (A, + A 2 ) 2 } x = 0, 

A 2T + [A 2XX + {A, + A 2 ) 2 } x = 0. (40) 

This case is corresponding to the resonances are located at { — 1, 2, 3, 4, 4, 6}. 
P-integrable model 5. 

Ayr + (A 1XX + A\) x + 2A 2 A 1X = 0, 

A 2T + {A 2XX + A\) x + 2A 1 A 2X = 0. (41) 
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This case is corresponding to the resonances are located at {—1, 0, 2, 4, 6, 7}. 
P-integrable model 6. 

A 1T + A lxxx + (A 1 + A 2 )(3A 1 + A 2 ) x = 0, 

A 2T + A 2XXX + (A x + A 2 )(3A 2 + A x ) x = 0. (42) 

This case is corresponding to the resonances are located at {—1, 1, 2, 4, 6, 6}. 

4. Exact solutions 

In this section, we study some types of exact solutions for the general couple KdV 
system (|27j) and (|2*%|) and some special types of P-integrable models. 

4.1. Travelling periodic and solitary wave solutions of the general coupled KdV 
system (l2~7jl and (|2~Hj) . In [6 , it is pointed out that some special types of exact solu- 
tions including travelling wave solutions of various nonlinear systems can be obtained via 
the deformation and mapping approach from the solutions of the cubic nonlinear Klein- 
Gordon equation (or namely 4 model). It is quite easy to see that for some types of 
travelling wave solutions of the coupled KdV system ()27|) and ()28|) can also be obtained 
by some suitable deformation relations from the travelling wave solution of the 4 model. 
For the travelling wave solution of the coupled KdV system (|27j) and (J2HJ), 

A, = A x (0 = A 1 {kX - kcT), A 2 = A 2 (0, (43) 

we have 

aiv4 1€ A 2 + {a 2 A\ + a z A x A 2 + a A k 2 A^ + a h A\ - cAi) € = 0, (44) 

5 X A^A 2 + {5 2 A{ + + ^k 2 A 2 ^ + 5 5 A 2 2 - cA 2 )^ = 0, (45) 

To map the travelling waves of the cubic nonlinear Klein-Gordon equation to those of 
the coupled KdV system, one may use different mapping relations such as the polynomial 
forms 6\, rational forms [7 j and may be more complicated and/or derive dependent forms 
jS]. However, here we just give a simple polynomial deformation relation 

A 1 =a + ai0(£) + a 2 (f)^) 2 , A 2 = b + M(0 + 6a 2 0(O 2 , (46) 
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where </>(£) is a travelling wave solution of the cubic nonlinear Klein-Gordon equation, 
i.e., satisfies 

4 = ^ + ^A0 4 + a (47) 

It is not very difficult to find that {A\, A 2 } given by ()46j) with (}4Tj) is a solution 
of the coupled KdV system (J2*7|) and (j2HI) if and only if the eleven solution parameters 
a , ai, a 2 , 6 , £>, A, C, fc, c and ten model parameters ««, i = 1, 2, 8 
satisfy the following eight constraints 

(2a! + 3a 3 + <oa 2 b)a 2 bi + [(30:36 + 6a 5 + a\b)a 2 + 3/c 2 a4A]a! = 0, 
{2a 3 b + 4a 5 )a 2 a + 2a 2 5 2 + («i + 2a 3 )ai&i + 2a 5 a 2 

+ [8£; 2 a 4 /i - 2c + (2a 3 + 4a 2 6 + 2«i)6 ]a 2 = 0, 
(2a 5 ai + b 1 a 3 )a + 2a 2 b bi + [/c 2 a 4 /i — c + («i + a 3 )6 ] a i = 0, 
(6i5 3 + 25 2 ai)a + {k 2 5^ + 260^5 - c)&i + (5i + 5 3 )6 «i = 0, 
(4<5 5 6 2 + 25 x b + 4<5 2 + 4<5 3 6)a 2 + 12fc 2 5 4 6A = 0, 
(4<5 2 + 2S 3 b)a 2 a + 2<5 5 6 2 + (2<5 3 + S^aih + 2b 2 a\ + [8/c 2 <5 4 6/i 

-2c6 + (A5 5 b + 25 x + 25 3 )6 ]a 2 = 0, (48) 
(Aa 3 b + 2ai6 + 4a 2 6 2 + 4a 5 )a 2 + 12fc 2 a 4 A = 0, 

[(65 5 b + 25! + 35 3 )a 2 + 3£; 2 5 4 A]&i + (35 3 b + 5 x b + 65 2 )a 2ai = 0. (49) 

It is obvious that the algebraic equation system ()49|) possesses many kinds of solutions. 

Here we just write down a most important and simplest solution when <5 4 = a 4 : 

7 7 n m2 6fc 2 Aa 4 . . 

a = a x = b = bi = 0, c = 4fc //a 4 , a 2 = — — — — — , (50) 

2a 5 + 2oa 3 + ba± + 2b z a 2 

while b is linked the model parameters by a cubic algebraic equation 

S 2 + ^5 3 - o 5 + l^ij b + (s s - a 3 - ^aij 6 - o 2 6 3 = 0. (51) 

More concretely, if we take </>(£) as the Jacobi elliptic conoid function 

(f) = cn(f , m) 
11 



which is a special solution of the (p A model with the parameters 



H = 2m 2 - 1, A = -2m 2 , C = 1 - m 



then we have a simple periodic wave solution for the coupled KdV equation (|27)l and (|2*5|) 
with 5 4 = «4 

where 6 is a solution of (|5ip. Furthermore, when m = 1, the periodic solution (|52|) 
becomes a simple solitary wave solution 

A x = 12fc2 " 4 — sech 2 (kX - Ak 3 a 4 T) , 

2a 5 + 2ba 3 + bay + 2b 2 a 2 V ' 

A 2 = 12^a4& sech 2 f kx _ Ak s T \ (53) 

2a 5 + 26a 3 + 6«i + 2b 2 a 2 V ; 

4.2. r-function solutions and Multi-soliton solutions of the coupled KdV sys- 
tem. 



4.2.1. The first type of r-function and multi-soliton solutions related to the KdV reduc- 
tions. It is straightforward to verify that for the coupled KdV equation system ()27|) and 
with 

5y = ay, (54) 

one can find at least one type of multiple soliton solutions because there is a simple KdV 
reduction 

A 1T + ayAyxxx + (aay + 2a 2 a 2 + 2aa 3 + 2a^)AyA lx = 0, A 2 = a Ay, (55) 
where a is a solution of the algebraic cubic equation 

2o?a 2 + {ay + 2a 3 - 25 5 )a 2 + (2a 5 -8y- 25 3 )a - 25 2 = 0. (56) 
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Then the coupled KdV equation system (J27)l and (J25j) with (J53J) possesses the following 
r function and multiple soliton solutions 

Ai = — = (lnr) vv , 57 

a aaci + 2a 2 a 2 + 2aa 3 + 2a 5 xx 

where r is just the usual r function. For the multi-soliton solution, the r function reads 

N N N 

T 1 ^ ^ TTj "I - ^ ^ ■Ai 1 i 2 7Ti 1 7li 2 -\- ^ ^ A-i-^i^iz ^11^12 ^13 ••■ -^1112. ..ijv ^"11^*12 •••""ijv ' 
1=1 «1<«2 «1<«2<«3 

7Tj = exp(/cjX - a 4 kfT), A hi2 ___ ik = j j A iaib . (58) 

j a <«i>,a,£>=l,2,...fc 

It is interesting that there is only one parameter condition (154)1 to get multiple soliton 
solution (}o7)|) while it has been proved that the model is non-Painleve integrable. In other 
words, the existence condition for multiple soliton solutions is not a sufficient condition 
of the integrability. 

Especially, because of there are three real solutions of (|5T)J) for the special coupled KdV 
equation (j37J), we can obtain three types of multiple soliton solutions {ui, Vi}, {u 2 , v 2 } 
and {u 3 , v 3 }, 

12 

v 1 = u 1 = — — (lnr)xx, (59) 

(ci - l)(c 2 - 1) 

12 

u 2 = t 7T7 -(lnr)xx, v 2 = CiU 2 , (60) 

(ci - l)(Ci - c 2 ) 

and 

12 

M 3 = 7 7T7 r(hir)xx, ^3 = c 2 u 2 (61) 

(c 2 - l){ci - c 2 ) 

with r being given by (|55j) . 

4.2.2. T/ie second type of r- function and multi-soliton solutions of the coupled KdV sys- 
tem. The multi-solitons of the coupled KdV system listed in the last subsection are ob- 
tained from its special KdV reduction. In [5| it has been found that even if for the 
non-integrable cases, the coupled nonlinear system may have much more abundant soli- 
tary wave structure. So we believe that for the coupled KdV system there may be other 
types of multiple soliton solutions. 
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For instance, if the model parameters have the following conditions 

aiO!2(«i<53 — 5ia 3 ) 7^ 0, <5 4 = 0:4, 

1 a 2 5i aiS 3 <5i«3 1 «i5 2 , . 

5s = 2 a3 + ^T- W a5 = _ ^7 + 2 53 + ^T' (62) 



"3 = T2 x — ( 63 ) 



25 2 aj _ + 63) 25 x a 2 
then we have a new type of multiple soliton solution 



12«ia 4 , a i >i / ft> i\ 

where r is still the r function of the usual KdV equation and for the multi-soliton solution 
it is given by (|58|) while A 2 is linked to the r-function by a linear equation 

^2r + 7 7 A 2X {hxr)xx + 7T "FT2 [( vxx b 

dia 3 - aid 3 (<)ia 3 - aid) 3 ) 2 

12aia 4 (5i5 3 + <5 2 + 25 2 a!) 

A 2 {\nT)xxx + OiAA 2 xxx = 0, (65) 



Sia 3 - aiS 3 

If the third condition ()63|) is not satisfied, then a nonlinear term 

r , r , <*i«3 , 25 2 ai 25fa 2 \ 
0i + d 3 -\ 1 1 A 2 A 2X 

Cti Oi Oil J 

has to be added to the left hand side of (|65|). 

Similarly, under the conditions 

zn r n a x "2(^3 -C1) 2 1 . . 

a 3 f: 0, «i = di = 0, d 4 = a 4 , d 5 = 7-73 -, (66) 



0=1(20^3 + 2<y 2 - eg) 

" 2 = 7^2 67 

2ci(ci - 5 3 ) 2 



where the constant C\ is determined by 



Cl = a 5 ± - 25 2 , (68) 
we have the following new type of multiple soliton solutions 

Ax = i lnr xx + 7-^^, 69 

Ci 5 3 - Ci 

where r is same as given by (|58*|) while v4 2 is also linked to the r function by a linear 
equation 

12a 4 (ci<53 + 2<5 2 ) r , „ N , 1445 2 (5 3 — ci)a 4 r/ , . . . 

A 2T + a 4 A 2XXX + ^[A 2 (lnr) xx ] x + 21 2 1J 4 [(lnr)^ x ] x = 0.(70) 
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In the same way, if the parameter condition (|B7jl is not satisfied, then we have to add a 
nonlinear term 

2a 2 (<5 3 -c 1 ) a 3 (2c l 5 3 + 25 2 -c 2 1 )\ 

1 77 s ^2^2X 

a 3 ci(5 3 -ci) / 

to the left hand side of (|70)l. 

4.2.3. T/ie £/iird type 0/ ' r -function and multi-soliton solutions of the coupled KdV system. 
Actually, in additional to the above types of multiple soliton solutions there may be 
other types of soliton solutions. Here is a further simple example for the following more 
specifical model 

A 1T + aA lxxx + bA x A xx + bcA 2 A 2X = 0, 

A 2T + aA 2XXX + bA x A 2X + bA 2 A lx = 0. (71) 

For this special model, its first type of multiple soliton solutions has the form 

A l =^0nr)xx, A 2 = ±±=A 1 (72) 
and the second type of multiple soliton solutions is given by 

Al = l^(lnr)xx ± V^A 2 , (73) 

while A 2 satisfies 

A 2T + aA 2XXX + 12a[A 2 (\n t) xx } x ± 2b^A 2 A 2X = 0. (74) 

For the special coupled KdV system (fTTj) . there is the following third type of multiple 
soliton solutions, 

A 1 = ^[Hr 2 1 +r 2 2 )] xx , (75) 

12a / t 2 \ 
A 2 = ± — -= arctan — , (76) 



xx 



where 



r = Ti + ir 2 

is just the usual r function of the KdV equation but with complex parameters, T\ and r 2 
are the real and imaginary parts of r respectively. 
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Figure 1. The interaction plots of two soliton solution for the field Al 
A 1 expressed by (J7SJ) and (|77jl at times (a) T = —10, (b) T = —5, (c) 
T = —0, (d) T = 5 and (e) T = 10 respectively. 

Fig. 1 and Fig. 2 are the special interaction plots of two soliton solution for the coupled 
KdV system (|7Tj) of the field A\ and A 2 expressed by ()75|) and (|75jl respectively with 

r = 1 + (1 + z)e felX - fc ? T + (1 + 3i)e k2X - k * T + {h - h) * Ui - 2 ) e (^)*-(*?+*f)r , 

(fci + fc 2 ) 2 

fci = l, = - (77) 
at times T = —10, —5, 0, 5 and 10 respectively. 
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Figure 2. The interaction plots of the two-soliton solution for the field 
A2 = A 2 expressed by (JZSJ) and (77|) at times (a) T = -10, (b) T = -5, (c) 
T = —0, (d) T = 5 and (e) T = 10 respectively. 



5. Summary and discussions 

In summary, a general type of coupled KdV system is derived from the coupled Euler 
equation system (JTJ) and (j2J) with (/5 ^ 0) and/or without (/3 ^ 0) the consideration of 
the earth rotation effects. In the derivation procedure, the frequently used inconsistent 
^/-average trick in the past literature is removed. 
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The integrability of the derived KdV system is checked by means of the well known 
Weiss- Tabor-Carnevale's Painleve test procedure. It is found that there are six types of 
Painleve integrable subcases for the derived coupled KdV system. 

The deformation and mapping method are used to get some types of travelling wave 
solutions including the conoidal periodic waves and single solitary waves for the general 
derived coupled KdV system with (|54p. 

It is found that the coupled nonlinear systems may possesses much more abundant 
solution structures. This phenomena is found before for the coupled non-integrable high- 
dimensional Klein-Gordon equation [£L In this paper we found that whence some kinds 
of model parameter conditions are satisfied, then there may be some different kinds of 
multiple soliton solutions and r function solutions. On the other hand, for the coupled 
KdV equation obtained here 

The dynamics of atmospheric blockings has been one of the central and important 
problem since they are the main representations of the general circulation anormaly in 
the areas of mid-high latitudes. Atmospheric blocking events have a strong influence on 
regional weather and climate. The observations have shown atmospheric blockings may 
locate in the mid-high latitudes, usually over the ocean, as the dipole pattern which was 
first discovered by Rex [TT] . 

For one layer atmospheric model, the single soliton solution of the constant coefficient 
KdV equation is responsible for the dipole pattern of the atmospheric blockings. To 
explain the blocking life cycle, one has to use the soliton solutions of the variable coefficient 
KdV equation ^2|- Using the similar analysis as the single layer atmospheric model, the 
soliton solution of the coupled KdV equation can also be used to explain the blockings 
under the two layer atmospheric description. In the same way, the soliton solutions of the 
constant coefficient coupled KdV equation proposed here can not be used to describe the 
blocking life cycle. To describe the blocking life cycle one has to extended the coupled 
KdV system (JTJ and (j2j) to variable coefficient case and the problem will be detailed 
studied in the near future study. 
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